Population Migration

100m National Shp

library(tidyverse)
Warning: package 'ggplot2' was built under R version 4.3.1
── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
✔ dplyr     1.1.2     ✔ readr     2.1.4
✔ forcats   1.0.0     ✔ stringr   1.5.0
✔ ggplot2   3.4.4     ✔ tibble    3.2.1
✔ lubridate 1.9.2     ✔ tidyr     1.3.0
✔ purrr     1.0.2     
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(sf)
Linking to GEOS 3.11.0, GDAL 3.5.3, PROJ 9.1.0; sf_use_s2() is TRUE
library(mapview)
The legacy packages maptools, rgdal, and rgeos, underpinning this package
will retire shortly. Please refer to R-spatial evolution reports on
https://r-spatial.org/r/2023/05/15/evolution4.html for details.
This package is now running under evolution status 0 
unzip("PWD_100m_national_SHP.zip")

pwd00 <- read_sf("PWD_100m_national_SHP/PWD_2000_national_100m.shp")
pwd05 <- read_sf("PWD_100m_national_SHP/PWD_2005_national_100m.shp")
pwd10 <- read_sf("PWD_100m_national_SHP/PWD_2010_national_100m.shp")
pwd15 <- read_sf("PWD_100m_national_SHP/PWD_2015_national_100m.shp")
pwd20 <- read_sf("PWD_100m_national_SHP/PWD_2020_national_100m.shp")
mapview(pwd00, col.region = "grey40") +
  mapview(pwd05, col.region = "#3099D3") +
  mapview(pwd10, col.region = "#924D7D") +
  mapview(pwd15, col.region = "#AB3967") +
  mapview(pwd20, col.region = "#DC143C")